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Abstract 

We use the Nikiforov-Uvarov method to calculate the bound states (energy spectra and wave 
functions) of a two-dimensional (2D) electron gas interacted with an exactly solvable pseudohar- 
monic confinement potential in a strong uniform magentic field inside dot and Aharonov-Bohm flux 
field inside a pseudodot. We give a unified treatment for both Schrodinger and spin-0 Klein-Gordon 
energy spectrum and wave functions as functions of chemical potential parameter, magnetic field 
strength, AB flux field and magnetic quantum number. We obtain analytic expression for the 
light interband absorption coefficient and threshold frequency of absorption as functions of applied 
magnetic field and geometrical size of quantum pseudodot. The temperature dependence energy 
levels for GaAs are also calculated. 
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I. INTRODUCTION 



Over a long time, a considerable interest has been paid for studying size effects in orbital 
magnetism [f ,2] and the magnetic properties of low-dimensional metallic and semiconduct- 
ing structures with restricted geometries [3] on nanostructures such as dots, wires, wells, 
antidots, well wires and antiwells [4,5,6]. Such structures can confine charge carriers in one, 
two and three dimensions. Experimental research is currently made to study the optical and 
quantum properties of low-dimensional semiconducting structures for the fabrication pur- 
poses and subsequent working of electronic and optical devices. More studies analyzing these 
structures have been focused on the interband light absorption coefficient in the spherical 
[7,8,9], parabolic, cylindrical and rectangular [10] quantum dots in the presence and absence 
of magnetic field [11]. More other works on optical properties in nanostructures [12,13], 
band structure calculations, transport properties of Aharonov-Bohm (AB) type oscillations 
[14] and Altshuler-Aharonov-Spivak (AAS) type oscillation [15]. 

The quantum antidot structure has been modeled in the presence and absence of repul- 
sive antidot potential, harmonic confining oscillator potential, the presence and absence of 
magnetic and Aharonov-Bohm (AB) flux fields in cylindrical coordinates [16]. This allows 
one to obtain an exact bound state solutions for the Schrodinger equation. The influence 
of dots and antidots on thermodynamic properties (e.g., magnetization) of the system, the 
magnetotransport properties and also the magneto-optical (MO) spectroscopic character- 
istics of a two-dimensional (2D) electron gas in a magnetic field are studied in [16]. The 
nature of MO transitions in this system demonstrate the appearance of rich spectrum of 
nonequidistant frequencies are different from the MO spectrum for a dot modeled by a har- 
monic confining potential. The quantum antidot is modeled as an electron moving outside 
a cylinder of radius a in the presence of magnetic and AB flux fields to find analytic expres- 
sions for energy and wave function [17]. The numerical and analytical solutions obtained 
for the dynamics of two classical electrons interacting via a Coulomb field in a 2D antidot 
superlatice potential in the presence of crossed electric and magnetic fields are quite dif- 
ferent than the noninteracting electrons [18]. Some authors have studied a 2D theoretical 
model for the quantum dot in which electrons were confined by a nonhomogenous magnetic 
field (the so-called magnetic antidot) [19]. The pseudoharmonic (PH) potential [20,21] is 
used in modeling the quantum dots (QDs) and quantum antidots (QADs) in nanostructures 
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[22]. The spectral properties in a 2D electron confined by a pseudoharmonic quantum dot 
(PHQD) potential in the presence of external strong uniform magnetic field il along the z 
direction in the presence of AB flux field created by a selenoid inserted inside the pseudodot 
have been studied. The Schrodinger and spinless Klein-Gordon equation are solved exactly 
for their bound states (energy spectrum and wave function) [22]. The advantage of the 
Klein-Gordon solution is that it provides us relativistic corrections to the commonly known 
nonrelativistic solution. 

It is well-known that factors such as impurity, electric and magnetic fields, pressure, 
and temperature play important roles in the electronic, optical and transport properties 
of low-dimensional semiconductor nanostructures [4,23-28]. In this regard, we carry out 
detailed exact analytic analysis of one-particle energetic spectrum and wave functions of both 
Schrodinger and Klein-Gordon equations with a pseudoharmonic potential in the presence of 
magnetic field and Aharonov-Bohm flux field by using the Nikiforov-Uvarov method [29,30]. 
The resulting energy spectrum serves as a base for calculating the corresponding interband 
light (optical) absorption coefficient and the threshold frequency value of absorption for 
the given model. In addition, the effect of the temperature on the effective mass is also 
calculated. 

The structure of the paper is as follows. In Sec. 2, the basic formulas of the Nikiforov- 
Uvarov (NU) method are outlined in short. In Sec. 3, we studied the nonrelativistic quan- 
tum dot and antidot with the pseudoharmonic potential in the presence of magnetic and 
Aharonov-Bohm flux fields. The exact analytic expressions for the energy spectra and wave 
functions are calculated. In Sec. 4, the analytic expressions for the bound states of the KG 
electron interacted via the pseudoharmonic potential in the presence of magnetic field and 
AB flux field are calculated. These basic formulas are also reduced to Schrodinger solutions 
for the pseudoharmonic potential model and free-field interactions under the non-relativistic 
limits. Results and discussions are performed in Sec. 5. The conclusions and outlook are 
presented in Sec. 6. 

II. NIKIFOROV-UVAROV METHOD 

This method is usually used in solving a second-order hypergeometric-type differential 
equations satisfying special orthogonal functions [29]. In spherical or cylindrical coordinates, 
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the resulting Schrodinger-like equation with a given potential is reduced to a hyp ergeo metric 
type equation through making a convenient change of variables, say, r — >■ s and then solved 
systematically for its exact or approximate eigensolutions (energy levels and wave functions). 
The most convenient equation, we consider here, takes the standard form 



where a(s) and a(s) are polynomials at most of second order, and r(s) is a first-degree 
polynomial and f(s) is a hypergeometric type polynomial. 

Next, we try to reduce Eq. (1) to a more comprehensible form by taking f(s) = (p(s)y(s) 
and choosing an appropriate function (f>(s): 

/(s) + um + m As) + (m + mm + m) y(s) = . (2) 

which appears to be more complicated than the standard form given in (1). To simplify (2), 
at first, we take the coefficient of y'(s), 

4>{s)^a{s) a(sY [) 



and set 

d/(s) tt(s) 



<f>(s) <t(s)' 

to obtain 



(4) 



*(s) = \[t(8)-t(s)], (5) 

where n(s) is a polynomial of degree at most one. Overmore, the above equation can be 
rewritten in the form: 

t(s) =t(s) + 2tt(s), (6) 

in which r(s) is a polynomial of order one. On the other hand, we can express the term 
<f>" (s) / <f>(s) appearing as one of the coefficients of Eq. (2) as 

4>"{s) _ (4>'{s)\ , {4/( 8 ) 



/ Tt(s)\' ( 7rfs) X 



<f)(s) V0( s )/ \(f>(s)J W s )/ \ a ( s ) 

In this case, the coefficient of y(s) can be simply recasted in the form: 



(7) 



<j)(s) cr(s) a 2 (s) a 2 (s) 



(8) 
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where 

a(s) = a(s) + tt 2 (s) + 7r(s)[f(s) - cr'(s)] + 7r'(s)<7(s). (9) 
Substituting the right-hand sides of Eq. (3) and Eq. (8) into Eq. (2), we finally obtain 

y'\s) + ^\y'(s) + ^-y(s) = 0. (10) 
a{s) a z {s) 

The above transformation allows one to set the hypergeometric function f(s) = (p(s)y(s), 
where 0(s) needs to satisfy the relation (4) with an arbitrary linear polynomial 7r(s). Thus, 
making the substitution: 

<t(s) = Xcr(s), 

where A is a constant. Hence, Eq. (10) turns into the so-called hypergeometric type equation: 

a(s)y" + T(s)y' + Xy = 0, (11) 

whose solution is already been given in [31]. Now, comparing Eq. (9) with Eq. (11) leads 
to the following quadratic equation: 

tt 2 (s) + [r(s) - a'(s)]n(s) + a(s) - ka(s) = 0, (12) 

where 

fc = A-7r'(s). (13) 
Thus, the solution of quadratic equation (12) is given by 



o-'(s)-t(s) , / / cr'(s) - t(s)\ 2 . , . . 

<s)= ' ±\ [ o -a(s) + ka(s), (14) 



where the parameter k inside the square root sign must be found explicitly to enable one 
to find the physical solutions of Eq. (14) for the plus and minus signs. Therefore, the 
expression under the square root sign has to be the square of a polynomial, since 7r(s) is a 
polynomial of degree at most one which provides an equation of the quadratic form available 
for the constant k. Having set the discriminant of this quadratic equal to zero, the constant 
k is determined clearly. Once the constant k is found, the task of the determination of the 
polynomial ir(s) from (14) becomes simple and straightforward. Further, r(s) and A can 
also be found from Eq. (6) and Eq. (13), respectively. 
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To make the solutions of Eq. (11) more general, we try to show that all the derivatives of 
hypergeometric type function are also of hypergeometric type. This can be easily acheived 
by differentiating Eq. (11) and letting i>i(s) = y'(s) 

a(s)v'l(s) + ri(s)«i(s) + p 1 v 1 (s) = 0, (15) 

where 7~i(s) = r(s) + cr'(s) and p 1 = X + t'( s )- r i( s ) * s a polynomial of degree at most 
one and p 1 is independent of the variable s. Equation (15) is obviously a hypergeometric 
type equation again. Further, taking v 2 (s) = y"{s) as a new representation and making the 
differentiation for the second time, we obtain 

a(s)v' 2 \s) + r 2 {s)v' 2 {s) + p 2 v 2 (s) = 0, (16) 

where 

r 2 (s) = n (s) + (/(s) = r(s) + 2a'(s), (17) 

A* 2 = t*i + = * + 2r '( s ) + ^"(«)- ( 18 ) 
Repeating this process, a general equation of hypergeometric type for v n (s) = y^ n \s) is 
constructed as a family of particular solutions corresponding to a given A; 

a(s)v'Xs) + r n (s)v' n (s) + pi n v n (s) = 0, (19) 

and hence the general recurrence relations for r„(s) and /i n can be found as 

r n (s) = r(s)+na'(s), (20) 

^ n = X + nr\s) + ^ Y ^a"(s), (21) 
respectively. When we set fx n = 0, then Eq. (21) becomes 

X = X nr = -nr'(s)-^ Y ^-a"(s), n = 0,l,2,... (22) 

and hence Eq. (19) has a particular solution 

v(°)=Vn(s) = -^±;[o n M8)], (23) 

which is known as the Rodrigues relation of degree n and p(s) is the weight function satisfying 

[a(r)p(r)]' = r{r)p{r). (24) 

Finally, to obtain an eigenvalue solution through the NU method, the relationship between 
X and X rir must be set up by means of Eq. (13) and Eq. (22). 
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III. NONRELATIVISTIC QDS AND QADS INFLUENCED BY MAGNETIC 
AND AB FLUX FIELDS 



A. Exactly solvable bound states 

Consider a two-dimensional (2D) single charged electron, e with an electronic effective 
mass, \i interacting via a radially symmetrical dot (electron) and antidot (hole). We will 
study the spectral properties of such dot and an antidot in a uniform magnetic field, ~B* = Bz 
and an AB flux field, applied simultanously. The Schrodinger equation is given by [32] 



, , , ■ - , . i>(?A) = Ei>(r,<l>), (25) 

where ~~X is the vector potential and the repulsive pseudoharmonic confinement quantum dot 
(PHQD) potential, V con f(r), describing the harmonic quantum dot and antidot structures, 
Vo(r) = Vor 2 jr\ and Vad(j') = Vor 2 /r 2 , respectively, is taken as [20,21] 

Konf (r) = V (- - , (26) 
\r r J 

where r and Vq are the zero point (effective radius) and the chemical potential. The 
vector potential A may be represented as a sum of two terms, such that 

^ x = if and ^ x ~~X 2 = 0, where if = Bz is the applied magnetic field, and ~X 2 
describes the additional magnetic flux <&ab created by a selenoid inserted inside the antidot 
(pseudodot). Hence, the vector potentials have azimuthal components given by [22] 

2 r ' 2nr \ 2 2vrr J V K J 

Let us consider the 2D cylindrical form of the wave functions: 

^(f,0) = -Le"^(r), m = 0,±l,±2,..., (28) 
v 2ir 



where m is the magnetic quantum number. Now, inserting the wave functions (28) into the 
Schrodinger equation (25), we obtain the following equation for the radial wave function 
g(r): 

9"{r) + V(r) + - ^ - 7 2 r 2 ) g(r) = 0, (29) 
where we have defined the parameters: 

^ 2 = f (E + 2V )-^(m + 0, (30a) 



(3 2 = (m + tf + a 2 , (30b) 

7 =*;? + br) ' (30c) 

where £ = §ab/&o with the flux quantum $ = hc/e, u c = eB / jic is the cyclotron frequency 
and a = kpr with hp = a/ 2fiVo / K 2 is the fermi wave vector of the electron. The magnetic 
quantum number m relates to the quantum number (3 [Eq. (30b)]. [1] Consequently, the 
radial wave function g(r) is required to satisfy the boundary conditions, i.e., g(0) = and 
g( r — > oo) = 0. In order to solve Eq. (29) by NU method, it is necessary to introduce the 
following variable s = r 2 , r G (0, oo) — t-sG (0, oo) which recasts Eq. (29) in the form of 
hypergeometric type differential equation (1) as 

9"(s) + ^g'(s) + (- 7 V + v 2 s - /3 2 ) g(s) = 0, (31) 

where we set g(r) = g(s). Applying the basic ideas of Ref. [30], by comparing Eq. (31) with 
Eq. (1) gives us the following polynomials: 

t(s) = 2, a(s) = 2s, a(s) = - 7 2 s 2 + v 2 s - (3 2 . (32) 

In the present case, if we substitute the polynomials given by Eq. (32) into Eq. (14), the 
following equality for the polynomial tt(s) can be obtained 

vr(s) = ±\Ji 2 s 2 + {2k-v 2 )s + (3 2 . (33) 

The expression under the square root of the above equation must be the square of a poly- 
nomial of first degree. This is possible only if its discriminant is zero and the constant 
parameter k can be determined from the condition that the expression under the square 
root has a double zero. Hence, k is obtained as k + _ = u 2 /2 ± j3^y. In that case, it can be 
written in the four possible forms of n(s); 

[ +{ 1S ±(5),ioTk + = \v 2 + M, 
tt(s) = { (34) 

{-( ri s±(3), for A;_ = \v 2 - p r 

One of the four possible forms of 7r(s) must be chosen to obtain an energy spectrum formula. 
Therefore, the most suitable form can be established by the choice: 



[1] For this system, only two independent integer quantum numbers are required. 
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for k-. The trick in this selection is to find the negative derivative of r(s) given in Eq. (6). 
Hence, r(s) and r'{s) are obtained as 

r(s) = 2(1 + /?)- 2 7 s, t'(s) = -2 7 < . (35) 

In this case, a new eigenvalue equation becomes 

A n = 2 7 n, n = 0,l,2,... (36) 

where it is beneficial to invite the quantity A n = — nr'(s) — ra ^~ 1 V // (s) in Eq. (22) with n 
is the radial quantum number. Another eigenvalue equation is obtained from the equality 
A = k_ + tt' in Eq. (13), 

A = y-7(/? + l). (37) 

In order to find an eigenvalue equation, the right-hand sides of Eq. (36) and Eq. (37) must 
be compared with each other, i.e., A n = A. In this case the result obtained will depend on 
E n ^ m in the closed form: 

v 2 = 2 (2ra + 1 + 13) 7. (38) 

Upon the substitution of the terms of right-hand sides of Eqs. (30a)- (30c) into Eq. (38), 
we can immediately obtain the following expression for the energy spectrum formula in the 
presence of PH potential : 

E n , m (£, 0) = hn(^n + ^f^J + \nw c (m + - 2V , Q = + (39) 

where \0\ = \J{m + £) 2 + a 2 > is an integer and cud = ^^Vo/fxr 2 . We have two sets of 
quantum numbers (n, to, (3) and (n', to', j3') for dot (electron) and antidot (hole), respectively. 
Therefore, expression (39) for the energy levels of the electron (hole) may be readily used 
for a study of the thermodynamic properties of quantum structures with dot and antidot in 
the presence and absence of magnetic field. 

If we ignore the last —2Vq term, the above formula becomes the Bogachek-Landman 
[16] energy levels, E n ^ m {^, (3) — Kl (n+ + ( m + > i n the presence of dot and 

antidot potential. In the absence of pseudoharmonic quantum dot (PHQD), i.e., Vq = 0, 
Q — > uj c , then E njTn (£) = ftw c [n + \{\m + £| + 1)] + \tko c (to + ^) which is the formula in 
the presence of magnetic and AB flux fields [16]. If we put £ = 0, i.e., in the absence 
of AB flux field, we find E n ^ m = hu c [n + |(|to| + to + 1)] which is the Landau energy 
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levels [33]. In the absence of magnetic field (uu c = 0) and an AB flux field (£ = 0), we 
find E n>m = (4hV /nr$) n+ (^m 2 + 2fiV r 2 /h 2 + l) /2 - 2V . When m = 0, we have 
E n = (AKVq/ hTq) (n + 1/2) for harmonic oscillator energy spectrum. 

Next, we calculate the corresponding wave functions for the present PH potential model. 
We find the first part of the wave function through Eq. (4), i.e., 

m ( S ) = exp(| ^ds^j =^I/V^ 2 . (40) 

Then, the weight function defined by Eq. (24) as 

p( s ) = _L e xp ( / 44^1 = s l/3| e" 7S , (41) 



which gives the second part of the wave function (Rodrigues formula); namely, Eq.(23): 

y n , m {s) ~ s ~ m ^J^ r (s n+m e-" s ) - ( 1S ) , (42) 

where (x) = { -^F (a,b + l ; x) is the associated Laguarre polynomial and Fa, b; x) is 
the confluent hypergeometric function. Using g(s) = m (s)y njm (s), in this way we may write 
the radial wave function in the following fashion 

9(r) = C ntm r^e-^ 2 I 2 F (-n, \fi\ + 1; 7 r 2 ) , (43) 

and finally the total wave function (28) becomes 



= .L , h m+i (" + ¥E r w\ e --<>*/i F (_„, m + 1; 7r 2j e , m * (44) 

\p\\ V 7in\ 

The energy levels in Eq. (39) differ from the usual Landau levels in cylindrical coordinate 
system [34] to which it transforms when £ = (i.e., <3>ab = 0), and a — > (i.e., when the 
chemical potential of dot and antidot vanishes, i.e., V — > 0). Nevertheless, the Landau 
levels are nearly continuous discrete spectrum for a particle confined to a large box with 
B — to equally spaced levels corresponding to B > 0. Each increment of energy, hu c , 
corresponding to free particle states, which is the degeneracy of each Landau level leading 
to a larger spacing as magnetic field B tends to become stronger [33]. The present model 
removes this degeneracy with energy levels spectrum becomes 



E n . m — t\U( 



1 .. 

n + - (\m\ + m + l) 



(45) 
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and the wave function reads as 



ml 



yn+1 ( n + m )! 



to! 



r m e - lr */2 F m + 1 . e im^ 



where 7 = (fiu c )/2h. In the limit when u c <^ g = w^-^then we have 



E nm — £0 + £iW c + £2^1 — £ 4^c + 



(46) 



(47) 



where 



e = -2V +N nm g, e x 



km 



-, £2 



N„ 



2g 



-, £4 



, TO + 1 ' 

« I n-\ — I , g 



B. Interband light absorption coefficient 



(48) 



Expressions (39) and (44), obtained above for charge carriers (electron or hole) energy 
spectrum and the corresponding wave function in quantum pseudodot under the influence 
of external magnetic field and AB flux field, allow to calculate the direct interband light 
absorption coefficient K(u) in such system and the threshold frequency of absorption. The 
light absorption coefficient can be expressed as [11,12,35]: 

K(p) = N E / #U,/^ m h n ,, m <A^ <^^# 6 (A - E^ - E*,^) , 

n m. ft m l Rl 



n,m,P n ' ,m' ,j3' 



n,m,P n',m',l3' 

( n + \P\)\(n' + 



T 



7r 2 n!n'!(|/3|!) 2 (|/3'|!) 2 



2vr 



e 

o 

/ 2\ |2 



. dre -(7+7')r 2 /2 r M + \P'\ 



x F (-n, |/3| + 1; 7 r 2 ) F (-n', + 1; 7V 2 ) | (A - E^ - E*^) , (49) 

where A = fuj — e g , e g is the width of forbidden energy gap, 00 is the frequency of incident 
light, N is a quantity proportional to the square of dipole moment matrix element modulus, 
ifj e(h) is the wave function of the electron (hole) and E e ^ is the corresponding energy of the 
electron (hole). 

Now, we use the integrals [33] 

2tt 



i(m+m' 



2tt if m = —m', 
if m 7^ —m', 



(50) 
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and 



J e~ KX x x ^F (-n, A; qx) F (-n', A; q'x) dx = T(\)K n+n '~ x (k - q)' 71 (k - q')~ n ' 
o 



x 2 Fi U,n',A 



(51) 



where is the Euler-Gamma function and 2 Fi (a, 6, c; z) is the hypergeometric function, 
to calculate the light absorption coefficient: 



n,m,/9 n',m',l3' 



where 



and 



Vn.n' 



aw 



(77' 



aI/8|+i 



7-7' 



rain'! 



/ 2 \ / 



477' 



(7 - y) 2 

Using Eqs. (39) and (49), we find the threshold frequency value of absorption as 



(52) 



(53) 



(54) 



(m + $W*o) + ^V rl/h 2 + 1 \ l/ gB V 8 y g/LB / $ AB 

= + o h/(-J +^+2^1™+^ 



/i'C / /iVg 2/i'c 



(55) 



When n = m = 0, then 
^00 = ^ s + 
ft 



+ 



- 4H- (56) 



C. Temperature dependence of the effective mass 

The variation of the effective mass with temperature is determined according to the 
expression [28,36,37] 



/'CO /('/') 1 + J ^ 



2 1 



(57) 
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where fj, e is the electronic mass, = 7.51 eV is the energy related to the momentum matrix 
element, A = 0.341 eV is the spin-orbit splitting and Eg(T) is the temperature-dependence 
of the energy gap (in eV units) at the V point which is given by [12,36,38,39] 

g r (r) = !.»!„_ (5-405 x 10-')^ 
3 v ; T + 204 v ; v ; 

Table 1 lists the temperature-dependent effective mass to the effective mass of donor electron, 
i.e., fi(T)/fi e for different values of temperatures. It is seen from Table 1 that raising the 
temperature will decrease the value of f(T) = fj,(T)/fi e . As a matter of fact, the decrease 
in this value means that kinetic energy of the donor electron decrease and consequently 
lowering the binding energy. The results are similar to Ref. [28]. Hence the temperature 
dependence energy spectrum formula can be expressed as 



En,m(B,T) 



. to 



2 



/CO 

which for GaAs turns to be 
E n ,m(B, T) = 14.9254^ 



J( m + £) 2 + a 2 /(T) + l\ 



l + 4-^/(T) \n+^ I 



-2V , 



I J(m + +0.067a 2 + l\ m + ^ 



1 + 0.268^ n + -1 + 



-2Vb, 



where we have used ji = 0.067 ji e . 

IV. THE SPINLESS KLEIN-GORDON PARTICLE IN MAGNETIC AND AB 
FLUX FIELDS 

The Klein-Gordon (KG) equation is wave equation mostly used in describing particle 
dynamics in relativistic quantum mechanics. Nonetheless, physically this equation describes 
a scalar particle (spin 0). Moreover, this wave equation, for free particles, is constructed 
using two objects: the four- vector linear momentum operator = ihd^ and the scalar 
rest mass M, allows one to introduce naturally two types of potential coupling. One is the 
gauge-invariant coupling to the four- vector potential {A^(~?)} 3 ^ =0 which is introduced via 
the minimal substitution P M — > — gA^, where g is a real coupling parameter. The other, 
is an additional coupling to the space-time scalar potential S COI1 f(l^) which is introduced 
by the substitution M — > M + S C(mi (~?). The term "four- vector" and "scalar" refers to 
the corresponding unitary irreducible representation of the Poincare space-time symmetry 
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group (the group of rotations and translations in (3 + l)-dimensional Minkowski space-time). 
Gauge invariance of the vector coupling allows for the freedom to fix the gauge (eliminating 
the non physical gauge modes) without altering the physical content of the problem. Many 
choose to simplify the solution of the problem by taking the space component of the vector 
potential to vanish (i.e., A). One may write the time-component of the four- vector potential 
as gA = V con i{r), then it ends up with two independent potential functions in the KG 
equation. These are the "vector" potential V con f(~?) and the "scalar" potential S con f(~?) 
[40,41]. 

The free KG equation is written as 

(9^ + M 2 )^ G (t,^) = 0. (59) 

Moreover, the vector and scalar couplings mentioned above introduce potential interactions 
by mapping the free KG equation as 

| c 2 + e _p^ + [M + S conf (r-)] 2 } tf?, </>) = [E- UOf 1>P>, 0), (60) 

where ^(~r^,0) is 2D cylindrical wave function defined as in (28). This type of coupling 
attracted a lot of attention in the literature due to the resulting simplification in the solution 
of the relativistic problem. The scalar-like potential coupling is added to the scalar mass 
so that in case when S con f(f) = ±V^ on f( r ), the KG equation could always be reduced to a 
Schrodinger-type second order differential equation as follows 



(f + -J) + 2(E±Mc 2 )V coni (l>) + M 2 c 4 -E 2 



^(T^,0) = O. (61) 



Hence, the bound state solutions of the above two cases are to be treated separetely as 
follows. 



A. The SconfC^) = +Vc 0nf (^) case 



The positive energy states (corresponding to S con f(r) = +V COQ {(7') in the nonrelativistic 



limit (taking E - Mc 2 = E and E + Mc 2 = 2/ic 2 , where \E\ < Mc 2 ) are solutions of 



ip(l>,<f)) = 0. 



(62) 



where ip(~?, 0) stands for either ^ + \~?, 0) or ip^ KG \l^, 0). This is the Schrodinger equation 
for the potential 2V con{ (l^). Thus, the choice S con f(~T^) = +V con{ (l^) produces a nontrivial 
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nonrelativistic limit with a potential function 2V con f(~?), and not K; 0n f(~r^). Accordingly, it 
would be natural to scale the potential term in Eq. (61) and Eq. (62) so that in the non- 
relativistic limit the interaction potential becomes Konf, not 2V con f. thus, we need to recast 
Eq. (61) and Eq. (62) as [41] 

2 



C 2 (f + -1) + (E + MC 2 ) Konf(^) + MV - E 2 



^,0) = O. 



(63a) 



(63b) 



with VconfCF") and ~X are given in Eq. (26) and Eq. (27), respectively. To avoid repeati- 
tion in solving Eq. (63a), we follow the same steps of solution explained before by taking 
^nra^! 0) = g(r)e tm< ^ /y/2n to obtain an equation satisfying the radial part: of the wave 
function: 



where we have used 

i2 _ 



9"{s) + j-^g'(s) + (-bjs 2 + \ 2 s - a\) g(s) = 0, 



1 [E 2 + 2{E + Mc 2 )V -M 2 c i ]-^(m + O, 



h 2 c 



bl = 



\ 2h 



h 2 C 2 TQ 



(E + Mc 2 ) V . 



(64) 

(65a) 

(65b) 
(65c) 



The solution of Eq. (64) can be easily constructed on making the changes: v — \ Ai, — > ai, 
and 7 — >■ b\. Thus, the equation for the KG positive energy states can be easily found from 
Eq. (38) as 

A 2 = 2(2n + l + a 1 )b 1 , (66) 
and further inserting Eqs. (65a)-(65c), we finally obtain the transcendental energy formula 



h\ l + 2n+ \ m' 2 + 



7 'I 



2M 



7i 



2uj' 2 ! 
u 2 c + -^ 7l = — ( 72 + 2V ) 7i - hu c m', 



7i = 



'E + Mc 2 ) 



, 7 2 = E - Mc 2 , u\ 



D 



2V , 2 _ 2MV r 2 

(X — 



c- - ' y Mrf " h 2 

where m' = m + £ is a new quantum number. We may find solution to the above transcen- 
dental equation as E = E^g- In the nonrelativistic limit ( , y 1 — > 2M and 7 2 — > E), the above 
equation can be easily reduced to the simple energy spectrum formula given in Eq. (39). 



(67) 
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Overmore, under the above parameters mapping, the 2D KG wave function can be found 
directly from Eq. (44) as 



= \/ ^^ r ' ai ' e " 6ir2/2L i ai) (fcir 2 )e^. (68) 
B. The S* conf (^) = -V conf (!>) case 

In this case, we follow the same steps of solution in the previous subsection: 

g"(s) + -^g'(s) + (-b 2 2 s 2 + \js - a 2 2 ) g(s) = 0, (69) 

where we have used 

A? A 2 = ^ [E 2 + 2 (E - Mc 2 ) V - M 2 c 4 ] - ^ (m + , (70a) 

„2 
_0 

ft 2 c 2 



a 2 -> a 2 = (m + 2 + [E - Mc 2 ) V , (70b) 



Thus, the equation for the KG negative energy states can be readily found as 

A 2 . = 2 (2n + 1 + a 2 ) b 2 , (71) 
which provides the transcendental energy spectrum formula 



(2n + 1) he + \l h 2 c 2 (m + £f + r 2 V ol2 } , / M 2 u 2 c + -^7 2 



'0 



= (c 2 7l + 2V ) 7 2 - hcMuj c (m + , (72) 
and the corresponding 2D KG wave function is found as 



Z>2| + 1 I 



^nlirA) = \l - g - |a2|)! r |a2| e^- 2 / 2 ^)(6 2 r 2 )e^. (73) 



It should be noted that the negative energy states (corresponding to S con t( r ) = — V con t( r )) 
are free fields since under these conditions Eq. (61) reduces to 

[ a+ ... \ 2 

2/i 



^, m (^,0)=O. (74) 
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which is a simple free-interaction mode. For the free fields, Eq. (74), the set of parameters 
in Eqs. (70a)- (70c) reads 



A 2 = \/-^ — (m + 0, a 2 = m + £, h = — , (75) 



which lead to the energy spectrum formula 



£ ( } = (n + m + €+^J hou c , 



(76) 



and and wave function 



V. RESULTS AND DISCUSSION 

We solved exactly the Schrodinger and Klein-Gordon equations for an electron under the 
pseudoharmonic interaction consisting of quantum dot potential and antidot potential in 
the presence of a uniform strong magnetic field if along the z axis and AB flux field created 
by an infinitely long selenoid inserted inside the pseudodot. We obtained bound state solu- 
tions including the energy spectrum formula (39) and wave function (44) for a Schrodinger 
electron. Overmore, for the Klein-Gordon electron, the positive energy equation (67) and 
wave function (68) is found for S con {(r) = +V con {( r ) case. However, the negative energy 
equation (67) and wave function (68) are found for S con {(l^) = +V con f(l^) case. These two 
cases are reduced to the Schrodinger equation with a potential interaction K; 0n f(~r^) and 
free field interaction solutions, respectively. Now we study the effect of the pseudoharmonic 
potential, the presence and absence of magnetic field B, the presence and absence of AB 
flux density £ and the antidot potential on the energy levels (39). To see the dependence 
of the energy spectrum on the magnetic quantum number, m, we take the following val- 
ues: magnetic field B — (6 T) AB flux field £ = 8, chemical potential Vq = 0.68346 



(meV) and r = 8.958 x lO" 6 cm [22]. Thus, we obtained a = ^2{iV r%/h 2 = 11.997702, 



2lod = ^/SVo/i^rl = 0.3280381 u c and hw = 1.05243ftw c [34], the dependence of the energy 
spectrum, (39) on the n and m is given by 



E, 



n,m 



yj{m + 8) 2 + 12 2 + l\ i 



1.05243 n + + - (m + 8) - 1.9678584, for B = 6 T. (78) 

huj c \ 2 12 
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where m = 0,±1,±2,... and n = 0,1,2, For the lowest ground state (n = 0): 



E , m /hu c = 1.05243 ( J(m + 8) 2 + 12 2 + 1 1/2 + (m + 8) /2 - 1.9678584, for B = 6 T. 



Overmore, to show the effect of magnetic field B on the energy spectrum, we take val- 
ues for parameters £ = 8, Vq = 0.68459 meV and r = 8.958 x 10 -6 cm [22], where 



a = y/2fj,V Q rl/h 2 = 12.007617 and Au 2 D = 8V /fir 2 = 0.120039 x 10 24 {rad/s) 2 , the de- 
pendence of energy levels (39) on the magnetic field becomes 



In Figure 1, we plot the pseudodot energy levels in the absence (presence) of pseudodot 
potential (i.e., V — — > a — (V ^ — > a = 12)) and in the absence (presence) of AB flux 
field <&ab (i-e., £ = (£ = 8)) as a function of magnetic quantum number m for B = 6 T. 
As demonstrated in Figure 1, the Landau energy states [33] (i.e., Vo = — > a = 0, £ = 
and £ = 8) are degenerate states (see, long dashed and dotted solid curves) for negative 
values of m, however, the pseudodot potential removes this degeneracy (case when Vq ^ 
— >■ a — 12), (see, solid and dotted dashed curves). In the absence of pseudodot potential 
(a = 0) and presence of AB flux field (£ = 8), the degeneracy still exists (long dashed line). 
It is found that the energy levels of PHQD potential are approximately equal the Landau 
energy levels for large absolute m values. However, they are quite different for small absolute 
m values (— 12 ^ m ^ 13 when £ = and —20 ^ m -< 5 when £ = 8). It is also noted that 
as the quantum number n increases (n > 0), the curves are quite similar to Figure 1 but 
the energy levels are pushed up toward the positive energy for all values of m. In Figure 
2 to Figure 7, we plot the magnetic field dependence of the ground state energy £o,m(£, a) 
(in units of meV) in the presence and absence of pseudodot potential and AB flux field for 
several values of magnetic quantum numbers m = 27, 35, 1, 0, —24 and —16, respectively. It 
is shown in Figure 2 to Figure 7 that pseudodot energy increases with increasing magnetic 
field strength. Further, in the absence of pseudodot potential, magnetic field in the positive 
z direction removes the degeneracy for positive m values. In these Figures, the behavior of 
pseudodot energy as function of the magnetic field B is shown in the presence of pseudodot 
potential and AB flux field (solid curves), in the absence of pseudodot potential and presence 




E n>m (meV) = 0.1157705v / £ 2 + 3.8803305 n + 





(79) 
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of AB flux field (dotted curves) and the absence of pseudodot potential and AB flux field 
(dashed curves). 

To investigate the dependence of the energy levels on temperature, we take the values of 
parameters: B = 6 T, f = 8, V Q = 0.68346 (meV) and r = 8.958 x KT 6 cm [22]. Hence, the 
temperature dependence of the energy levels (in the units of Huu c ) at the T point are given 
by 



E n ,rn(T) 
fVjJc 



/CO 



\Jl + (0.32804) 2 /(T) n + 



m + 8) 2 + 144/(T) + 1 \ m + 8 



+ 



- 1.9678584, 



(80) 



where f(T) is calculated in Table 1 at any temperature value. In GaAs, we have f(T) = 0.067 
[11]. Taking the special values of parameters £ = 8, V = 0.68459 meV and r = 8.958 x 10~ 6 
cm [22], two parameters (temperature and magnetic field) dependence of the energy levels 
(in units of meV) are calculated as 



E n , m (B, T) 



/CO 



0.1 157705 ^B 2 + 3.8803305/(T) n + 



m + 8) 2 + 144/(T) + 1 



+ 0.11577055 



m + 8 



- 1.36918 (units meV). 



(81) 



which becomes 



E n>m (B) = 14.9254 



0.1157705V5 2 + 0.26 n + 



m + 8Y + 9.648 + 1 



+ 0.11577055 



m + 8 



- 1.36918 (units meV). 



(82) 



for GaAs. Figure 8 to Figure 12 show the variation of the pseudodot energy levels (in 
meV) as function of magnetic field B (in Tesla) (82) in the presence of pseudodot potential 
and AB flux field (solid curves), in the absence of pseudodot potential and presence of AB 
flux field (dotted curves) and the absence of pseudodot potential and AB flux field (dashed 
curves). for various values of radial quantum numbers n and magnetic quantum numbers m. 
For GaAs case, we consider the following cases (a) n = m = 0, (b) n = 5, m = 0, (c) n — 0, 
m = 5, (d) n = 0, m = — 5 and (e) n — 5, m — —5 in Figures 8 to 12, respectively. 



19 



VI. CONCLUSIONS AND OUTLOOK 



In this work, we have obtained bound state energies and wave functions of the KG particle 
in the field of pseudoharmonic quantum dot and antidot structure in the presence of a uni- 
form magnetic and AB flux fields. The positive (negative) KG energy states corresponding 
to S con f(l^) = +V con {(l^) (S con i(l^) = — V con {(~?)) are studied. Overmore, the Schrodinger 
bound state solutions are found. Under nonrelativistic limit, the KG equation with equal 
mixture of scalar and vector potentials S con {(~?) = +V con f(~?) and S con f(~?) = — V con f(l^) 
can be easily reduced into the well-known Schrodinger equation of a particle with an interac- 
tion potential field and a free field, respectively. Overmore, the nonrelativistic electron and 
hole energy spectra and the their corresponding wave functions are used to calculate the the 
interband light absorption coefficient and the the threshold frequency of absorption. Also, 
the energy spectrum of the electron (hole) may be used for a study of the thermodynamic 
properties of quantum structures with dot (antidot) in a magnetic field. The temperature 
dependence of the energy levels are calculated using the Table 1 at any temperature T 
(Kelvin). 
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TABLE I: Calculated f(T) with different values of temperature for GaAs. 



T (K) fi(T)/fi e T (K) M (T)/ Me 



n 


n nfifiQQS4 


1 7n 


fl 0653679 


i n 


n nfifiQssfi 


1 80 


D6591 77 




n nfifiQfins 

U.UUUc/UUO 


-L yu 


U.UUtJUU^O 




n OfifiQI 74 


909 


n n64QD8n 


4D 


n nfifisfins 


91 n 


fl 0647490 


uU 


n nfifi7Qi i 

u.uuu (yii 


99D 


D 084^874 

U.UU^iUO / 4: 


fin 

UU 


n nfifi7i 1 9 

U.UUU / lli 


9^0 


D 0fi449' : 5 c i 

U.UU^i^i^OtJ 


7D 


n nfififi9i 7 


949 


fl D649573 


80 




950 


0640891 a 


on 

cf\J 


n nfifi4i 78 


960 


D D6391 88 


1 nn 


n D663D51 


970 


fl D637468 


i in 


n nfifii 861 


98D 


fl D63573D 


120 


0.0660614 


290 


0.0633976 


130 


0.0659315 


300 


0.0632206 a 


140 


0.0657968 


350 


0.0623154 


150 


0.0656577 a 


400 


0.0613818 a 


160 


0.0655147 


500 


0.0594513 a 



a See Ref. [28]. 

FIG. 1: Pseudodot n = energy levels (in tko c unit) as a function of magnetic quantum number 
m in the presence and absence of PHQD potential (a = 12 and a = 0) and in the presence and 
absence of AB flux field (£ = 8 and £ = 0) for magnetic field B = 6 T. 
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FIG. 2: Ground state pseudodot energy levels (in meV) as a function of magnetic field B (in 
Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB 
flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence of 
AB flux field, respectively for magnetic quantum number m = 27. 



FIG. 3: Ground state pseudodot energy levels (in meV) as a function of magnetic field B (in 
Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB 
flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence of 
AB flux field, respectively for magnetic quantum number m = 35. 



FIG. 4: Ground state pseudodot energy levels (in meV) as a function of magnetic field B (in 
Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB 
flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence of 
AB flux field, respectively for magnetic quantum number m = 1. 



FIG. 5: Ground state pseudodot energy levels (in meV) as a function of magnetic field B (in 
Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB 
flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence of 
AB flux field, respectively for magnetic quantum number m = 0. 
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FIG. 6: Ground state pseudodot energy levels (in meV) as a function of magnetic field B (in 
Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB 
flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence of 
AB flux field, respectively for magnetic quantum number m = —24. 



FIG. 7: Ground state pseudodot energy levels (in meV) as a function of magnetic field B (in 
Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB 
flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence of 
AB flux field, respectively for magnetic quantum number m = — 16. 



FIG. 8: Ground state GaAs pseudodot energy levels (in meV) as a function of magnetic field B 
(in Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of 
AB flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence 
of AB flux field, respectively for magnetic quantum number m = 0. 



FIG. 9: GaAs pseudodot energy levels (in meV) as a function of magnetic field B (in Tesla). 
Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB flux field, 
Landau energy levels in presence of AB flux field and Landau levels in the absence of AB flux field, 
respectively for n = 5 and m = 0. 
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FIG. 10: Ground state GaAs pseudodot energy levels (in meV) as a function of magnetic field B 
(in Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of 
AB flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence 
of AB flux field, respectively for magnetic quantum number m = 5. 



FIG. 11: Ground state GaAs pseudodot energy levels (in meV) as a function of magnetic field B 
(in Tesla). Solid, dotted and dashed curves represent the pseudodot energy levels in presence of 
AB flux field, Landau energy levels in presence of AB flux field and Landau levels in the absence 
of AB flux field, respectively for magnetic quantum number m = — 5. 



FIG. 12: GaAs pseudodot energy levels (in meV) as a function of magnetic field B (in Tesla). 
Solid, dotted and dashed curves represent the pseudodot energy levels in presence of AB flux field, 
Landau energy levels in presence of AB flux field and Landau levels in the absence of AB flux field, 
respectively for n = 5 and m = — 5. 
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